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Abstract Here I discuss possible relations between 
free precession of neutron stars, Tkachenko waves inside 
them and glitches. I note that the proposed precession 
period of the isolated neutron star RX J0720. 4-3125 
(Haberl ct al. 2006) is consistent with the period of 
Tkachenko waves for the spin period 8.4 s. Based on 
a possible observation of a glitch in RX J0720. 4-3125 
(van Kerkwijk et al. 2007), I propose a simple model, in 
which long period precession is powered by Tkachenko 
waves generated by a glitch. The period of free preces- 
sion, determined by a NS oblateness, should be equal 
to the standing Tkachenko wave period for effective en- 
ergy transfer from the standing wave to the precession 
motion. A similar scenario can be applicable also in the 
case of the PSR B1828-11. 
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1 Introduction 

Isolated neutron stars (NSs) being non-spherical bodies 
are expected to demonst rate free pre cession (for a brief 
review see, for example, iLinkl 20031 ). However, exam- 
ples of this phenomena are less than few, and even in 
rare cases when a precession-like behavior is observed 
different interpretations can be dis cussed (even not re- 
lated to precession, see for example iRuderman and Gil 

The problem of long period free precession in NSs 
is a long standing one. A NS can precess if it is non- 
spherical and rotation axis does not coincide with a 
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principal axis. Typically, biaxial objects are discussed, 
so deviation from spherical symme try can be described 
by one parameter - oblateness (see lAkgun et al. ( 2006 ) 
for a discussion of triaxial model). Expected values of 
NS oblateness (due to rotation or influence of strong 
magnetic fields) can naturally lead to precession peri- 
ods about one year. The precession period is equal to 
-Pprcc = -P/e- Here P is the spin period of a NS, e - 
its oblateness, and P pre s _ precession period. Measured 
precession periods require oblateness about 10 -8 . 

However, discussing dynamics of NSs it is necessary 
to take into accout the network of supcrfluid vortices in- 
side them. The neutron supcrfluid liquid in the interior 
of a NS participates in rotation via formation of quan- 
tized vortex lines. The density of these lines per unit 
area is n = 2fl/k. Here = 2tt/P is spin frequency, 
and k = h/2m ri , where m n is a neutro n mass (see, for 
example. IShapiro and Teukolskvlll983l . Ch. 10). The 
vortices exist in the core of a NS, where they can in- 
teract with supcrfluid (superconducting) protons and 
normal electrons, and in the crust, where they can pin 
to it. 

Coupling of supcrfluid neutron vortices with clcc- 
trons in a core results in damping of free precession 
( Alpar and Ogelmanlll987l ). But the time scale of this 
damping is long enough, according to these authors. 
For spin period about 1 se cond it is ~ 400 - 10 4 P pre c 
( Alpar and Ogelmanl 1987 ). Still, this time scale is 
much shorter than a NS age, so some excitation mech- 
anism is necessary for precession. As it is discussed 
below, in the presented model excitation is due to a 
glitch. 

A kind of pinning ( "immobilization" ) of vortices can 
also happen in the core due to interactions with mag- 
netic flux tubes (see discussion, for example, in ILinkl 



2007). In this case, the moment of inertia of "pinned" 



neutrons (which is about I) is about 10 times larger, 
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than the moment of the remaining parts of a NS, I c . 2 Tkachenko waves 



So, P 



prcc 



0.1P. 



A different kind of problem appears if pinning in 
the crust is taken into account. For absolute pinning 
no long period precession is possible. Instead, the pe- 
riod of precession becomes equal to P(I/I p ), here / is 
NS moment of inertia, I p is the moment of inertia of 
pinned superfluid i n the crust. Ty pically it is expected 
that I p /I ~ 10~ 2 (|Shahamlll977r ). and the precession 
period is just ~ 100 P if the absolut e pinning is valid. 
However, lAlpar and Qgelmanl (|19871 ) showed that this 
is not the case due to finite temperature. Because of 
thermal effect always there is vortex creep which allow 
the pinned superfluid to follow precession. 

The best example of a NS with precession-like behav- 
ior is PSR B1828-11. The proposed period is about 51 1 
days with a harmonic at 256 days (|Stairs et al. 2000| ). 
Most of discussions related to free precession deal with 
this source. In particular, the problem of non-existence 
of long period precession for strong pinning is typically 
confronted with observaions of PSR B1828-11. 

Recently, appeared another possible example of long 
period free precession in NSs. The existence of ~ 7 
years precession period in one of a small group of 
isolated NSs (called XDINS - X-ray Dim Isolated 
NSs, or ICoNS - Isolated Cooling NSs, or the Mag- 
nificent Seven) - RX J0720. 4-3125 - was suspected 



(jHaberl et alJbooet Fl So, this object was added to the 
list, and the paradoxical situation of long precessi on in 
presence of s uperfluid vortices was reconsidered by [Link 
(200G. 120071 ). This author proposed that either protons 
in NS interiors are type I superconductors, or neutrons 
in the o uter core are normal ( i.e., n ot superfluid) . More 
recently Glampedakis et all ( 2008 ) demonstrated that 
for long spin period and small precession angles NSs can 
have long precession periods (note, that for PSR B1828- 
11 the preces sion angle is prop osed to be small, about 
few degrees (|Stairs et alj|2000j) . but for RX J0720.4 



3125 it can be larger > 10 degrees (IHaberl et aL. 
20061) ). So, a ccording tolGlampedakis et all |2008l ). the 



conclusion by Link! ( 20061 ) and other authors that in the 



strong drag regime P p 



0.1P can be under doubt 



due to a short wavelength instability. 

Clearly, the problem of free precession in NSs is far 
from being solved completely. In this brief note, based 
on coincidence between Tkachenko wave period and 
precession period in cases of PSR B1828-11 and RX 
J0720. 4-3125, 1 discuss a mechanism to support preces- 
sion in isolated NSs. 



A slightly different period ~ 4.3 years was proposed by van 
Kerkwijk and Kaplan (2007) based on timing analysis. However, 
these authors consider a model with a glitch to fit better due 
to a rapid change in spectral properties, see van Kerkwijk et al. 
(2007). 



A simple model for long period precession of isolated 
NSs pr oposed here is rel ated to the so-called Tkachenko 
waves |Tkachenkolll966l) . These are displacement waves 
in the vortex line array that exist in rotating superfluid, 
or in other words a kind of sound waves propagating in 
the lattice of neutron votices perpendicular to them. 
A good introduction to the Tkachenko waves physics 
can be found in the paper bv lAndereck and Glaberson 
( 19821 ). 

Already in early 70-s this phenomena wa s suggested 
to explain periodic modulations in NSs (| Ruder man 
197dt Dyson 1971 ). At that time motivation had 
been related to reported wobbling of the Crab pulsar, 
which was not confirmed by later observations. Then 
this approach was nearly forg otten, and only recently 
Noronha and Sedrakianl (|2008[ ) returned to considera- 
tion of Tkachenko waves in NSs. In particular, they 
demonstrated that behavior of PSR B1828-11 can be 
explained by the se waves. 

According to Rudermanl ( 1970t ) (see also an exam- 
ple given by Dvsonl ( 1971 )) the period of a standing 
Tkachenko wave in a NS can be estimated as: 



P T = (2w/k)(l/V T ) ~ 1.77P 6 P 1/2 yi's. 



(1) 



Here Vt - wave velocity, which in a simple case depends 
only on the spin period (Vt ~ P" 1 / 2 ) and fundamental 
constants. Rq - the core radius normalized to 10 6 cm 
(with such normalization the equation provides an es- 
timate close to the upper limit for the period). The 
estimate is made for the mode with kR = 5. Spin pe- 
riod of a NS, P, is given in seconds. 

As one can see, this period is of order of those related 
to free precession. 



3 Scenario for long period precession 

The proposed precession periods for RX J0720.4- 
3125 and PSR B1828-11 are very similar to the 
Tkachenko wave periods for these stars. For RX 
J0720. 4-3125 the peri od of precession is proposed to 
be equal to ~ 7 years (Haberl et al. 20061 ) or 4.3 years 
( van Kerkwijk and KaplanT 2007? ) . spin period of this 
NS is equal to ~ 8.4 s (see, for example, Haberll ( 2007 ) 
for a review on XDINS). For PSR B1828-11 the preces- 
sion period is equal to ~ 500 days, while the spin period 
is equal to 0.4 s, for this object the coincidence between 
precession and Tkachenko wav e periods was already 
mentioned (see, for example, iGusev and Kitiashvili 
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20041 iNoronha and Sedrakianll2008h . For R about few - 
10 km one obtains that the precession period is consis- 
tent with Pt for both NSs (also the mode can be used 
as a parameter, however everyw here here I use kR = 5, 
as proposed by lRudermarJ[l970T) . 

Moreover, in the case of RX J0720. 4-3125 pulse pro- 
file modulations and spectral changes are observed. Pe- 
riod modulations related just to the Tkachenko waves 
alone hardly can be responsible for such evolution. Pre- 
cession is necessary. The idea, proposed here is the 
following: energy stored in standing Tkachenko wave 
can power precession of a NS. The necessary condition 
for effective energy transfer can be the equality of Pt 
and free precession period. The latter one depends on 
oblatcncss of a NS, which can be due to strong mag- 
netic fields. The former one depends only on the spin 
frequency. Coincidence between the two characteristic 
time scales which depend on different quantities should 
be not a very frequent occasion. However, one more 
condition is necessary - it is necessary to generate the 
standing wave. The necessity of these two conditions 
ca n explain why pre cessing isolated NSs are so rare. 

Rudermanl (|l970h notes that Tkachenko waves can 
be generated by starquake glitches. Energy of pre- 



cession is E, 



prcc 



"r>r Bcfl 2 /2 , where 9 is the pre- 
cession amplitude (jjonesf \20oi ). For PSR B1828-11 
this energy is about 3 10 36 erg. This value is sig- 
nificantly smaller than typical glitch energy. In the 
case of RX J0720 . 4:-3125 a glitch was proposed by 
van Kerkwiik et al.l (|2007l ). The glitch energy was es- 
timated by them as ~ 10 37 5 erg (Av/v ~ 510 -8 ). 
Of course, it is necessary to say, that the energy of 
the glitch was estimated according to the fit with cu- 
bic model (jvan Kerkwiik et al.ll2007l ) . and so the value 
can be different if after the glitch the evolution is due 
to precession. Still, as an order of magni tude estimate 
the value from van Kerkwijk et alJ (J2007|) can be used 



Haberl et al. ( 2006| ) estimated the amplitude of pre- 
cession > 10°, but warn about uncertainties of their 
model. The period of precession is about 7 years. The 
energy of the precession motion then appears to be 
- 3 10 35 (6»/10°) 2 erg. With these values in hand the 
energy of the glitch is enough to drive precession even 
for large 9 if efficiency of energy transfer is not very 
low. 



van Kerkwijk et al.l |2007j) relate a "jump" in spec- 
tral properties of RX J0720.4-3125 to the glitch. In 
my opinion, this means that it is more probable that 
the glitch was due to a quake, not due to vortex lines 
unpinning (or accretion episode etc., see below). 

Taking altogether, for RX J0720.4-3125 the following 
scenario is proposed: a glitch (most probably due to a 
starquake) generates Tkachenko waves; the period of a 



standing Tkachenko wave is equal to the free precession 
period for this NS; due to the standing wave precession 
starts after a glitch, or just there is an energy input into 
the pre-existing precession motion. 

Tkachenko waves periodically change the spin fre- 
quency and moment of inertia of a NS. Waves move 
perpendicular to the vortex lines, which are parallel to 
the spin axis. The moment of inertia of a star can 
be non-symmetric respect to this axis, for example if 
oblateness is due to strong magnetic field. I speculate 
that periodic modulation of spin frequency and all com- 
ponents of moment of inertia in resonance with the pre- 
cession period (determined by oblateness) would lead to 
energy transfer from Tkachenko waves to the precession 
motion. 



4 Discussion 

Absence of free precession in absolute majority of iso- 
lated NSs indicates that this phenomena needs some 
rare coincidence in properties of a NS. Here it is pro- 
posed that it is necessary to have: 

• Pt ~ Pprcc, 

• a glitch to generate Tkachenko waves. 



Instead of the proposal by iLinkl |2007j) - "A slowly - 
precessing neutron star cannot glitch" - I propose 
another: slow precession is powered by glitches via 
Tkachenko waves. 

Observations of RX J0720. 4-3125 are roughly con- 
sistent with this scenario. On the other hand, in the 
case of PSR B1828-11 no glitches have been observed. 
However, it is necessary to study for how long pre- 
ces sion can survive aft e r a g litch. If an old estimate 
bv lAlpar and Ogelmanl |l987l ). 400 - 10 4 P proc is valid, 
then this time is long enough. If precession is peri- 
odically excited by glitches via Tkachenko waves even 



damp ing on a time scale of few precession cycles ([Link 



mp 
36) 



2006) would not contradict observations of RX J0720.4- 



3125. 



The glitch in RX J0720.4-3125 reported bv lvan Kerkwiik et al.| 
( 20071 ) by its consequencies is similar to the one ob- 
served in an ano malous X-ray pulsar (AXP) CXOU 
J164710. 2-455216 |lsrael et aljl2007l : lMuno et all2007h . 
After the glitch the luminosity of the source was in- 
creased, and its spectrum changed. So, the jump 
in properties of the spectr um and luminosity of RX 



J0720.4-3125 proposed by Ivan Kerkwiik et all (|2007l ) 



can be directly related to a glitch, which is weaker than 
in the case of CXOU J164710. 2-455216 (still similarities 
in behavior of these sources can be considered as a kind 
of support to the hypothesis of a link between AXPs 



•1 



and XDINS). But the evolution of the NS parameters 
after the "jump" requires precession. 

Note, that the timing solution before MJ D 52821, 

when a possible " glitch" happened according to van Kerkwijk et 
( 20071 ). can be relatively well described by the so-called t, 
cu bic solution (the second derivativ e of v is non-zero), 



one obtains that the reccurence time between two suc- 
cesive glitches is: 



sec 



van Kerkwijk and Kaplan! (|2007). Spectral changes 



before this date are not very la rge ( Haberl et al. 20061 ; 



van Kerkwijk and Kapl an 2007). After MJD 52821 the 



timing s olution is well described b y a pe riodic func 



van Kerkwijk and Kaplan ( 20071 ) (these au- 



tion, see 

thors studied several models with and without a glitch 
in their t wo papers), and sp ectral changes follow this 
law, too ( Haberl et al. 2006| ^1. Based on that, I sug- 
gest that the timing residuals might be also explained 
by a model without (or with small) precession before 
the glitch, and strong precession after. However, this 
particular model has never been tested quantitatively 
against observational data. 

Glitches naturall y can produce thermal afterglows 
(jHirano et all Il997h . About 10 38 - 10 43 ergs can be 
released in a glitch (in the case of RX J0720. 4-3125 
according to estimates of the inc rease in spin fre- 
quency bv Ivan Kerkwijk et al. 2007 this value is closer 
to 10 38 erg). However, Hirano et al. showed that a 
thermal response of a NS to a glitch cannot produce 
a smooth temperature increase on the time scale of 
years. If surface temperat ure is increased just by few 
percent, as it is required by van Kerkwijk et all ( 2007 ). 
then the brightening lasts just for few days (this corre- 
sponds to weak energy release). If we require a tem- 
peratur e rise for a long t ime, then the effect is too 
strong ( Hirano et al. Il997l ). So, I conclude that spec- 
tral changes on a long time scale should be attributed 
to precession of the NS. 

Glitches of AXPs (and soft gamma-repeaters) can be 
different in nature with respect to radio pulsar glitches, 
as the former can be related to crust fracture due to 
superstrong magnetic field. Still, the origin of a glitch 
is not important for our discussion here. "Normal" 
glitches are quite common for long period pulsars, for 
example, PSR J1814-1744 with spin period about 4 
secon ds demonstrated a glitch (jJanssen and Stappers 
20061 ). So, RX J0720.4-3125 can glitch not only via 



the mechanism operating in magnetars, but also due to 
convenient mechanisms proposed for normal radio pul- 
sars. Fo r them one can estimate the reccurence time 
following I Alpar and Bavkall (|l994l ). 

If the glitch of RX J0720.4- 3125 is due to unpinning , 
then using standard formulae (jAlpar and Bavkallll994l ) 

2 However, van Kerkwijk and Kaplan (2007) propose that the 
period is not close to ~7 years, but is ~4.3 years. Still, on a 
short time scale - since MJD 52821 - this is not very certain. 



(2) 



Parameter <50 is the critical value of the difference be- 
tween the rotation frequencies of normal matter and 
the superfluid at a boundary layer. <5fJ itself can be 
estimated as (Alpar & Baykal 1994): 



Afi— , 



(3) 



here J p is the effective moment of inertia of the region 
of a pinning layer. 

Combining these two formulae one obtains the rela- 
tion for the time between glitches: 



2 1 AD, 
in ~n 



■t oc t, 



(4) 



where t is the age of a pulsar, t = fi/2 Cl. 

Surprisingly, the time is about 10 years for RX 
J0720. 4-3125. I.e. it is quite probable to observe one 
since the discovery of this object. Then, one can expect 
to see similar phenomcnae in other Magnificent seven 
objects. However, they are less studied, and may be 
some glitches are missed. 

Still, to produce luminosity and spectral changes, 
and to generate Tkachcnko waves, it is probably more 

natural to have a glitch due to a starquak e. In the 

case of the quake model (jAlpar and Bavkall Il994l ) the 
time between glitches is longer, about 300 years for RX 
J0720.4-3125. 



2(A + B)(f)An/n 



(5) 



The est imate above was obtain ed assuming standard 
values (jAlpar and Bavkall Il994l) A = 10 52 erg, (f> = 
10" 3 , B = 10 48 erg, and I ~ 10 45 g cm 2 

Then, we can be just lucky to find a glitch in ~ 10 
years of observations (but note, that it is not the only 
XDINS observed). Or, in XDINS quakes do not follow 
th e formula for radio pulsars . 

van Kerkwijk et al.l ( 20071 ) proposed that an accre- 



tion episode can be responsible for spectral changes af- 
ter a "glitch" in RX J0720.4-3125. I think that this is 
not a very probable reason for the origin of the glitch 
and corresponding changes. It is hardly possible to 
imagine that if we observe such an episode just after 
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10 years of observations, other episodes were not fre- 
quent during the evolution of this source. With frequent 
episodes of accretion of light elements a NS should fol 



low a slightly different cooling track (jKaminker et al 



20061) . Such stars with accreted envelopes are hotter in 



their youth, but colder after they mature. This is more 
similar to the properties of 1E1207.4-5209 and Kes 79 
( Gotthelf and Halpernll2007l ). 

A remarkable difference between some NSs in su- 
pernova remnants (so-called CCOs) and XDINS in the 
solar vicinity can be due to the existence of accreted 
envelopes in the former. Note, that wc do not ob- 
serve descendants of 1E1207-Iike sources in our prox- 
imity. If they follow a standard cooling curve, like RX 
J0720.4-3125 and other XDINS, then they have to be 
observed. On the other hand, we do not see ancestors 
of XDINS in supernova remnants. This can be related 
to the fact that CCOs descendants are too cold at the 
age of XDINS to be easily detected by ROSAT, and 
vice versa ancestors of XDINS are not hot enough to 
be easily found in some supernova remnants. 

In this note I neglect (as most, if not all, other au- 
thors who studied Tkachenko waves in NSs) the influ- 
ence of interaction between neutron vortex lines and 
magnetic flux tubes. This interaction can significantly 
affect the velocity of waves, and so their period, and to 
damp themU The velocity of a Tkachenko wave can be 
estimated as: 
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Vt = l/2(MV27TTO pair ) 1/2 = (ftfc/87r) 1/2 ~ b/P. (6) 

Here TO pa i r = 2m n , and b is the distance between vortex 
lines. So, Pt ~ (R/b)P in the simple case when there 
are no interactions with flux tubes or other complica- 
tions. When vortices are able to "communicate" with 
the help of numerous magnetic flux tubes, the velocity 
of the wave can be larger, so the period of Tkachenko- 
like wave would be shorter. This question should be 
explored. 

To conclude, in this short note I proposed that long 
period precession in RX J0720. 4-3125 can be related 
to Tkachenko waves, generated in a recent glitch. In 
the proposed model before a glitch precession could be 
negligible. The critical condition for the free preces- 
sion excitation is the equality between Tkachenko wave 
period and the period of free precession. 
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